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Abstract. In a previous analysis, Maurin et al. (2001) have constrained several parameters of the cosmic ray diffusive 
propagation (the diffusion coefficient normalization Kq and its spectral index S, the halo half-thickness L, the Alfven 
velocity Va, and the galactic wind Vc) using stable nuclei. In a second paper (Donato et al., 2001), these parameters 
were shown to reproduce the observed antiproton spectrum with no further adjustment. In the present paper, we 
extend the analysis to the /3-radioactive nuclei ^°Be, ^^Al and ^^Cl. These species will be shown to be particularly 
sensitive to the properties of the local interstellar medium (lism) . As studies of the lism suggest that we live in an 
underdense bubble of extent rhoie ~ 50 — 200 pc, this local feature must be taken into account. We present a modified 
version of our diffusion model which describes the underdensity as a hole in the galactic disc; we believe some of the 
formulae presented here are new. It is found that the presence of the bubble leads to a decrease in the radioactive 
fluxes which can be approximated by a simple factor exp{—rhoie/lrad) where Irad = K^tq is the typical distance 
travelled by a radioactive nucleus before it decays. We find that each of the radioactive nuclei independently points 
■ towards a bubble of radius < 100 pc. If these nuclei are considered simultaneously, only models with a bubble radius 
^sj I Thole ~ 60 — 100 pc are marginally consistent with data. In particular, the standard case Vhoie = pc is disfavored. 
• Our main concern is about the consistency of the currently available data, especially ^^Al/^^Al. 
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1. Introduction 

The idea of using long-lived radioactive isotopes to estimate the average age of cosmic rays is not new and the name 
"radioactive chronometer" is now standard. Indeed, the finite lifetime of these nuclei has a deep impact on their 
propagation and they are expected to give independent clues about the propagation parameters. To this end, the 
nucleus ^°Be {ti/2 = 1.51 Myr) was proposed in the fifties, but due to an insufficient knowledge of the processes 
occurring during propagation, it has long been thought not to be the best candidate after all. Instead, Reames (1980) 
proposed ^'^Mn which is EC (electronic capture) unstable {ti/2 = 3.74 Myr), but it was subsequently shown to be of 
no great interest, owing to the specific behavior of the EC channel. Later on, Casse (1973) proposed ^"^Mn which is 
a mixed /3-EC unstable as a chronometer of the Fe group, the latter being mostly created by collisions of cosmic ray 
Fe on the interstellar medium. These propositions are very instructive since they demonstrate the early hesitations 
and progressive comprehension of cosmic ray nuclei propagation. From an observational point of view, the first clues 
about this propagation time came from the isotope ^''Be (Hagen et al., 1977; Webber et al., 1977; Garcia-Munoz et 
al., 1977). Other isotopically resolved chronometers followed as ^®A1 (Freier et al., 1980), "^^Cl (Wiedenbeck, 1985) and 



^^Mn (Leske, 1993). Since then, numerous accurate experiments have been performed (see section 5.2) giving more 
confident estimates of the surviving fraction of these radioactive clocks. 

Most of the conclusions of radioactive studies have been drawn in leaky box models, using the notion of surviving 
fraction to estimate the age of cosmic rays. The early analyses focused on the surviving fraction and this quantity hinted 
at the existence of a diffusive halo far beyond the thin galactic disc (see in particular Simpson and Garcia-Munoz, 1987 
for a discussion). We emphasize that in diffusion models, these notions have no clear physical meaning (see discussion 



in section 3.1), as stable species and each radioactive species have different ages. The main goal nowadays is rather 
to use a diffusion model that can consistently reproduce all observed radioactive abundances. Our purpose here - 
following Ptuskin et al. (1997), Ptuskin and Soutoul (1998) - is to demonstrate that one may have to consider local 
properties of the interstellar medium (ism) in order to bring this program to a successful conclusion. 

The plan of the paper is the following: in a first part we recall the nuclear properties of the radioactive species. In the 
second part, we focus on the astrophysical aspects of propagation. More specifically, we introduce an inhomogeneous 
model for the production and propagation of radioactive species. We also discuss the meaning of older models and 
some false claims that are made in the literature about the predictive power of /3 radioactive clocks. We then present 
our analysis and results. 



2. Description of radioactive species 

From the propagation point of view, one can associate three attributes to each nucleus: (i) a reaction cross section, 
(ii) a spallation production cross section and (iii) a rest frame half-life. We will not discuss further the first two 
points as they were discussed at length in Maurin et al. (2001) (hereafter Paper I). Indeed, the cross section code and 
parameterizations developed in Webber et al. (1990), Webber et al. (1998) (spallation cross section) and Tripathi et 
al. (1997a; 1997b; 1999) (reaction cross section) are appropriate for both stable and radioactive nuclei. 

However, it is not so straightforward to specify decay channels for unstable nuclei. Despite the fact that one can 
find half- lives easily in nuclear charts (Audi et al., 1997)^, the information needed from a cosmicist point of view is not 
always clearly extractable: for example (3 and EC modes are often combined whereas these two effects act differently 



^ see also 



http://sutekh.nd.rl.ac.uk/CoN /I or |http:/ /nucleardata.nuclear. lu.se/nucleardata/ 
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during propagation. Thus in the foUowing section, we give a hst of nuclei that need to be propagated as unstable. We 
emphasize that these ingredients were already present in the version of our code used in Paper I, although they were 
not explicitly discussed. 

2.1. Half -lives 

Two different disintegration modes must be distinguished: f3 decay (the nucleus decays spontaneously with a given 
lifetime tq = ii/2/ln2) and EC decay (the nucleus decays after capturing a K-shell electron if one has been attached). 
In the latter case, the decay rate depends on the electron density and on the typical attachment time. On Earth, 
many electrons are available so that a nucleus unstable via electronic capture behaves as Z?"*" (at least in terms of the 
daughter nucleus and half-life). Nevertheless in the iSM, where fewer electrons can be attached, the attachment time 
must be taken into account in the EC mode. We will come back to this point later on. 

Among all the radioactive nuclei, those which have lifetimes of the order of the propagation time may bring 
interesting information about their propagation. We select and classify them into three categories: /3-unstable, EC- 
unstable and mixed /3-EC-unstable. We emphasize that this classification only has a meaning in the context of cosmic 
ray propagation: were we interested in stellar nucleosynthesis, for example, other tables would be necessary to encode 
disintegration channels as the characteristic times of evolution involved and the properties of the medium would be 
drastically different. 

Up to now, the most complete description of half-lives adapted to cosmic ray description has been given in Letaw 
et al. (1984). In this paper, we only consider (3 unstable nuclei, but as results about EC species will be presented in a 
separate paper (Donato et al. 2001, in preparation) we find it convenient to describe in this section all unstable species 
(/3 and EC). Thus, the following tables are similar to those in Letaw et al. (1984), but we restrict ourselves to Z < 30 
and we take the opportunity to incorporate updated values for some half-lives. 

2.1.1. /3 unstable 

In Table H, we present the three pure /? unstable isotopes having proper half-lives in the time range 1 kyr-100 Myr. 
Nucleus C is the shortest-living /3 unstable isotope included in our cascade. This choice is motivated by the measured 
propagation time, which runs between 10-20 Myr depending on the propagation model (this time decreases with 
energy). If a nucleus has a half-life < kyr, we can consider it as a ghost nucleus (see details in Paper I); this means 
in particular that as soon as it is created, this nucleus immediately (compared to the typical propagation time) 
disappears into its daughter nucleus, so that its density in cosmic radiation is negligible and it does not need to be 
propagated. We might object our tq > 1 kpc cut-off, arguing that because of the Lorentz factor 7 = Etot/'m, even 

Table 1. Pure /3 unstable isotopes (1 kyr < ti/2 < 100 Myr). 



z 


Nucleus 


Daughter 


t5'72"- (error) 


4 
6 

26 


iOBe 
14 

6 

IgFet 


^^N 
(i?Co'^)iNi 


LSI"''" (0.06) 
S.TS*^" (0.04) 
(0.3) 



' The transition from Fe to Co has a half-life of 1.5 Myr while transition from Co to Ni is immediate from a cosmic ray point 
of view {ti/2 ~ 5 yr). 

a short proper lifetime can give a large effective lifetime at sufficiently high energy. Actually, we are concerned with 
energies Etot < 100 GeV/nuc, so that 7 < 100. It follows that proper half- lives tq < 1 kyr correspond to effective 
half-lives r — 7T0 < 0.1 Myr, which is too short compared to the propagation time for our purpose. Moreover, this 
cut-off is natural in practice as there are no radioactive half-lives in the interval 0.1-1 kyr. Conversely, if this lifetime is 
too high, the corresponding nucleus could be propagated as a stable one. To give an example, '*°K - whose f3 half-life 
is about ^ 10^ yr - can clearly be considered as stable. 

2.1.2. EC unstable 

In Table ||, we present the isotopes that are unstable under electronic capture transitions. Unlike for /9-decay, these 
nuclei have to attach an electron to be able to decay. As (i) nuclei are completely stripped of e~ at cosmic ray energies 
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and (ii) the interstellar medium is very poor in e~, the attachment time may be much longer than the lifetime of the 
decay. This means that even a nucleus with a very short EC-lifetime may have a long effective lifetime. As a result, 
there is no need to consider any lower bound on the half-lives. This will be discussed further in section 2.2, and the 
interested reader can refer to Letaw et al (1984) or Adams et al (1985) who show in a leaky box model that the 
effective half-life of EC nuclei is at least of the order of the attachment rate, which is about ~ Myr~^. 



Table 2. Pure EC unstable isotopes. 



z 


unstable (ec) 


Daughter 


tlY^ - (error) 


4 


iBe 




53.29'* (0.07) 


18 


?lAr 




35.04'' (0.04) 


20 


iCa 


Sk 




22 


ttTi 


(!tSc^)!^Ca 


49'"' (3) 


23 


iv 


iTi 


330'' (15) 


24 


tICrt 


(iv^-:)iiTi 


21.56'' (0.03)t 


24 


51ri-. 

24<or 




27.702'' (0.004) 


25 


iiMn 


53 

24 or 


3 74Myr 


26 


iFe 


iiMn 


2.73^" (0.03) 


27 


i?Co 


ia^Fe 


271.79'' (0.09) 


28 


28^^! 


iCo 


80'=^'' (11) 



^ This nucleus has an allowed ji transition, but contrary to ^'*Mn and ^''Ni (cf Table it has not been studied recently so that 
we can set it as a pure EC. 

* In this two-step reaction, the second transition 48y/^48rpj -^^^ a half-life greater than the first one - 15.9735'' (0.0025). 
Nevertheless, this second reaction can be taken as immediate because of its fi nature. We thus can consider this second element 
as a ghost (see text). Finally, only the first reaction (''^Cr— >■**¥) enters the decay rate. 

^ This nucleus has a fi decay, but with ti/2 > 100 Gyr. For the same reason as explained in section [2.1. l[ , it is sufficient to take 
into account the EC channel. 



2.1.3. Mixed /3-EC unstable 

Finally we turn to the case in which the two decay modes are allowed (hereafter mixed decay) , in table ^. In general, 
/3 decay is dominant, but we also need to consider some EC contributions. This is the case for ^"^Mn which is often 



Table 3. Mixed EC-/3 isotopes 



z 


Nucleus 


Daughter (ec) 


t\Y^- (error) 


Daughter [fi) 


tlY^- (error) 


13 
17 
25 
28 


?iAi* 

??C1 
ilMnt 
iNi* 


?iMg 

36 o 
16^5 

54 p 

24 or 
^(PCo^-)iFe 


4.08"'"' (0.15) 

312.3'' (0.4) 
6.075'' (0.020) 


?iMg 

36 Ar 

ilFe 
§(i?Co^)iFe 


0.91"'"' (0.04) 

Q 3Q7Afl/r (0.002) 

0.494"'"' (0.006) 
0.051'^^'"' (0.022) 



* Martmez-Pinedo and Vogel (1998) (for this nucleus, one can often find in the literature a half-life of 0.74 Myr: it corresponds 
to the combined /3-EC channels, and this is incorrect since these two processes are not equivalent during propagation). 
^ (3 half- life taken from Martmez-Pinedo and Vogel (1998). A higher value for this channel is not excluded, 
t Lund Fisker et al. (1999). 

^ 27C0 decays via EC (80%) and (20%), but as the half-life is of the order of two months, one can consider that the only 
effective channel is /3-decay so that this nucleus vanishes immediately. Notice that these values are taken from Goldman (1982). 
More recent references (Audi et al., 1997) or nuclear charts on the web (see footnote 1) are ignored because they give either 
pure /3 channel or pure EC channel. 
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discussed as a possible chronometer for the Fe group, but its EC contribution is rarely treated correctly. Incidentally, 
table ^ also shows that the usual two chronometers ^®A1 and '^^Cl have mixed decay channels. But contrary to ^^Mn 
we will see in the next section that the EC channel is completely suppressed, and it is usually neglected. 

2. 2. Modifications of decay properties during propagation: effective half-lives 

While (3 decay is independent of the nucleus environment, electronic capture decay depends sharply on the iSM 
electronic density as well as on the attachment and stripping cross sections. It also follows that the energy dependence 
of lifetimes is more complex in the EC mode than in the /3 mode. 

This is of great importance for the stability of nuclei, and it gives a flavor of how different the propagation is 
for these two modes, even for similar half-lives {i.e. tEc = t/3)- In order to clarify these specific behaviors, the pure 
/3-decay will serve as a reference process to which EC processes can be compared. The expression for the density of a 
/3-unstable nucleus can be found in Webber, Lee and Gupta (1992), or formulae (A8), (A9) and (AlO) of Paper I (the 
disintegration rate appears explicitly in Si). 

As regards the other radioactive species, they actually may be treated like /? radioactive species as it is possible to 
transform their solution into that of pure /3 with an effective rate instead of the pure EC rate or of the usual combined 
/3-EC rate. An effective lifetime can be introduced, obtained by combining tEc, attachment and stripping cross sections 
(and tp for mixed decay). A discussion of the validity of such approximations can be found in Letaw et al. (1984), 
Adams et al. (1985) and references therein. The physical inputs required to describe EC-unstable and /3-unstable are 

Table 4. Pure (3 unstable isotopes (1 kyr < ii/2 < 100 Myr) from a propagation point of view. 



z 


Nucleus 


Daughter 


tS'/a"- (error) 


4 

6 
13 
17 
26 


i^Be 

6 

iAl 
f?Cl 


i°B 

?iMg 

fiAr 


1.51"'"' (0.06) 
5.73''^'' (0.04) 
O.gi""" (0.04) 
Q 3Q^Mar (Q 002) 

1.5"'"' (0.3) 



different and we made the choice here to focus on the latter case. The case of EC-unstable species will be discussed 
elsewhere (Donato et al. 2001, in preparation). Consequently, the nuclei that one might consider as pure /? are given 
in Table |. 

We have checked that the EC mode can be neglected for ^^Al and ^^Cl but not for '^'^Mn and ^^Ni. The comparison 
of our tables to Letaw et al. (1984) shows that the largest difference is for ^"Fe half-life which changed from 3 x 10^ 
to 1.5 X 10^ yr. For other nuclei, minor corrections in half-lives and their uncertainties have been made. 

Finally, we recall that the propagation of nuclei which are radioactively produced requires a specific treatment. 
Their solution corresponds to formula (All) of Paper I and as they do not play any role in this study, we will not 
discuss these daughter nuclei further. 

3. Motivation for a local treatment of radioactive cosmic rays propagation 

The propagation of cosmic rays is a priori affected by the details of the gas density distribution in the galactic disc. 
Several models have been proposed to take this into account: for example. Strong and Moskalenko (1998) consider a 
gas density distribution n{r,z), whereas Osborne and Ptuskin (1987) and Ptuskin et Soutoul (1990) model the ism 
as a cloudy medium. As far as the stable species are concerned, one can use an equivalent treatment with an average 
description of the ism density. This is what we did in Paper I. 

The situation is drastically different for radioactive species. First, due to their finite lifetime, those that reach the 
solar system must have been created locally, i.e. in a region which is at most a few hundred parsecs away (see below). 
As a consequence, the local interstellar matter (lism) has to be carefully described. 

Actually, there are reasons to believe that this lism is highly inhomogeneous, which motivates a more elaborate 
model for radioactive species. This model will be described in the next section, after we have discussed the locality of 
radioactive production and the inhomogeneity of the lism. 
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3.1. Locality from diffusion equation 

It is often claimed that the size of the diffusion halo can be estimated from radioactive cosmic ray species. This is 
true in the framework of a leaky box propagation model, but not for diffusion models (see below). Indeed, the leaky 
box models have been so widely used in the past that some of their conclusions have become popular wisdom and 
arc sometimes used out of the proper context. Actually, as was explicitly shown in the early seventies (Prishchcp and 
Ptuskin, 1975; Ginzburg et al., 1980), leaky box models are almost never equivalent to diffusion models, so that one 
should be very careful before directly applying leaky box inspired results to another class of models, e.g. diffusion 
models. 

3.1.1. The old leaky box paradigm, and what we should forget to go further. . . 

The so-called leaky box modelling of cosmic ray propagation assumes that the particles freely move in a homogeneous 
finite-sized box. When a particle reaches a boundary of the box, there is a finite probability that it escapes the system. 
As a result, forgetting about anything but propagation, the density is given by, in this model. 



where Tcsc{E) is the typical time a particle of energy E spends in the box. It is possible to account for the measured 
abundance of all the stable cosmic ray nuclei abundances with a suitable Tcsc{E), adjusted for all nuclei. However, the 
physical picture provided by this model is wrong. Cosmic rays do not freely stream in a homogeneous box, but they are 
scattered by the inhomogeneities of the galactic magnetic field. This results in a diffusive propagation. This difference 
is of no importance when one is concerned only with the local abundance of stable nuclei, and it can be shown that 
in this case, there is an equivalence between these two approaches (see discussion in Paper I and references therein). 
But as soon as one tries to change the mathematically correct local description of stable nuclei provided by the leaky 
box into a physical picture and tries to extend it either to radioactive species or to another location, trouble begins. 

First, it is obvious that from the start, the leaky box model is unable to describe the spatial distribution of cosmic 
rays. We will not discuss this point further. 

Second, it is also clear that this model is bound to fail for radioactive cosmic ray nuclei. The species we consider 
have a proper lifetime of the order of Trad ~ 1 Myr, so that at the speed of light they can travel more than ~ 300 kpc 
before they decay. As the radius of the galaxy is ~ 20 kpc, these nuclei should have enough time to propagate in 
the whole available volume and their distribution is sensitive to the global geometry and localization of the sources. 
This is at variance with the diffusion model in which the average distance travelled by a nucleus during a time Trad is 
given by ^,/K{E)^\^, which is of the order of a few hundred parsecs at low energy. This means that in this case, the 
distribution of these species only depends on the local characteristics of the diffusive medium and of the sources. 

It is sometimes claimed that a radioactive species gives the size L of the diffusive halo: the bigger the halo, the 
smaller the quantity of nuclei that survive from the sources to the Earth. This is easily understood in the framework 
of the leaky box, but it is physically wrong for diffusion models. A more quantitative version of the intuitive argument 
given above can be found in Prishchep & Ptuskin (1975). They showed that diffusion models are equivalent to leaky 
box models only when the relationship tq » L'^/K is satisfied. In the most favorable of the cases we consider here, 
this corresponds to tq 3> 100 Myr, which is wrong for the nuclei we are interested in. 

Thus, one should be careful not to use leaky box induced conclusions out of this specific context. They may be 
very intuitive but wrong when applied to a diffusion model. For example, the next section shows that the radioactive- 
to-parent ratio is not sensitive to the size of the halo in a diffusion model. 

3.1.2. Does a radioactive really tell something about the halo size L ? 

In the diffusion model, a radioactive nucleus with a proper lifetime tq has most probably travelled a distance Irad = 
\J K{E)-jTo between its creation in the disc and its detection on Earth. At low energy, this distance is much smaller 
than the size of the diffusive volume, and this nucleus is not very sensitive to the boundaries of this volume. This can 
be seen in the expression giving the ratio radioactive/parent nucleus. In the case when there is only one parent nucleus 
and if energetic gains and losses are discarded, the density is given by (see equation A8 of Paper I) 



dN{E,t) _ N{E,t) 

dt Tesc{E) 



= 




(1) 



with 



Ai = 2/ir"'=' +Vc + KSi coth(5,L/2) 



(2) 
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and 



1/2 



(3) 



^^2 + ^2 + j 

When the lifetime is so short that T^L/K L/lrad ^ 1, we have coth(5'iL/2) « 1. It can be easily seen that in 
these conditions, the ratio N'^^''{z = 0, r)/A^P'"'(z = 0,7") becomes independent of the halo size L. Figure |l| displays 
this ratio for various L when formula (|^) is used. The diffusion coefficient K[E) and l^ad are increasing functions of 
energy, so that when Irad becomes of the order of the typical size of the halo {L or R), diffusive propagation is affected 
by the shape and size of the halo. 



N""=(Rs„,.z = 0)/N^°-"XRs„„z = 0) 



0.03 



.025 - 



0.02 



.015 - 



0.005 



L=100 kpc 



''" L=1 kpc 



L = 3 kpc 



L= 1 kpc 



Energy [GeV/ n uc] 



Fig. 1. unstable/parent ratio ( °A1/ °Si) for different values of L in kpc, all other parameters being fixed. 



We have presented here a toy model in which only one parent was considered. Actually, there are several parents 
and the above conclusion is not exact but remains approximate if they are considered, so that at low energy, the ratio 
j>^uns I ^^^^ ]\[pa-r jg |-j^g same whatever the size of the halo. All current data on radioactive nuclei are at low energy, 
so that we do not expect such a ratio to give direct constraints on L. 

Nevertheless, the ratio of isotopes having the same parents, like ^"Be/^Be, are usually considered. As the quantity 



10 



^par ^^"^^ does depend on the halo size the ratio °Be/^Be eventually also does 



3.2. Local sub-density: we all leave in a local bubble 

Interest in the study of our local environment has grown in the last thirty years (see for example Cox et Reynolds 
(1987) for a review). As emphasized in McKee (1998), the very recent all sky survey by rosat became an invaluable 
resource for all astronomers in general, and for the lism study in particular. In this section, we present evidence 
indicating why it is probably necessary to have a specific description for the local iSM (in a region of about ^ 100 pc 
around the solar neighborhood) . We then develop a model to incorporate the effect of these local properties of the iSM 
on radioactive species within our diffusion model. 
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3.2.1. Main properties of the lism 

First, it has to be noticed that far more progress has been made in mapping three dimensional distribution of galaxies 
many megaparsecs away than the distribution of local interstellar clouds within a hundred parsecs (McKee, 1998). 
We can nevertheless summarize a few points about the local composition. First, the lism is defined as a region of 
extremely hot gas (~ 10^ — 10^ K) and low density {n < 0.005 cm~^) within a bubble of radius between < 65-250 
pc surrounded by a dense neutral gas boundary ('hydrogen wall') (Sfeir et al., 1999; Linsky et al., 2000). A smaller 
scale description of this bubble shows that the Sun is located in a local fluff with Nhi ~ 0.1 cm^^, T ~ 10"* K and a 
typical extension of 50 pc. It is of great importance for further modelling to realize that the local bubble is highly 
asymmetric (Cox et Reynolds, 1987; Welsh et al., 1994; Fruscione et al., 1994), and that several cloudlets are present 
in the bubble (for a schematic representation, see for example fig. 1. of Breitschwerdt et al., 2000). Various models 
have attempted to explain the formation of this local bubble (Smith and Cox, 2001; Breitschwerdt et al., 2000; see 
also Cox, 1997 for a brief review), but this subject is far beyond our concern. 

Finally, a model can be built, considering that the Sun is surrounded by a first shell of radius ^ 50 pc and density 
~ 0.1 cm^'^, and a second shell of radius ~ 200 pc with an almost null density. Beyond this second shell, we recover 
the usual average density 1 cm~'^ (or zero density if the radius of the second shell extends beyond the disc). In the rest 
of this paper, we use a coarser description, in which the underdensity is modelled with one cylindrical hole of radius 
''hole to be determined and with a null density. It will be called 'the hole' throughout this paper. 



3.2.2. Consequences for the radioactive production 

The propagation of cosmic rays is influenced by (i) the local magnetic properties and (ii) the local matter content of 
the disc. In this section, we investigate the effect of this hole on these properties. 

First, we assume that diffusion itself, as described by the coefficient K{E), is not affected by the presence of the 
hole, i.e. diffusion is homogeneous. Homogeneity seems to be confirmed by radio and 7-ray observations, which can 
test in situ the spectrum and density of cosmic rays (McKee 1998, MorfiU & Freyberg 1998). 

The presence of the hole has another consequence. Because the density is lower, there are less spallations in the 
bubble, so that the secondary nuclei abundances (including radioactive ones) are probably perturbed. As a matter of 
fact, a realistic description of the matter content of the disc should take into account a spatial distribution (exponential 
or more complex) as used for example in Strong et Moskalenko (1998). However, the only relevant quantity for stable 
nuclei is the average grammage of matter they cross. As these nuclei propagate in a region which is much bigger than 
the hole, this average grammage is only slightly changed by the presence of the hole. This can be seen in figure ^ 
where the change in the radial distribution of a stable species due to the hole is represented. In the following, the local 
density of stable nuclei will be computed with a full matter disc of density 1 cm~'^. 

The situation is different for radioactive species, as the typical distance they travel from their creation in the disc 
to their detection on Earth is limited by their finite lifetime r. This distance can be estimated as Irad ~ y/ K{E)'^to. 
For the nuclei we consider here, and for energies of a few hundred MeV/nuc, Irad is of the same order of magnitude 
as the size of the hole, so that their propagation is expected to be more strongly perturbed. To be specific, the fact 
that spallations do not occur within the hole has three distinct physical effects. First, it leads to a decrease in the 
spallation source term of the radioactive species. Second, it also leads to a local decrease of destructive spallations. 
Third, as there is less interstellar matter to interact with, the energy losses are also lowered. 



4. Modelling the local propagation 

We now want to take these remarks and incorporate the three physical effects discussed above in our diffusion model. 
We start with a very simple model, in which the influence of the hole can be easily understood. We then adapt the 
model discussed at length in Paper I and Donato et al. (2001) (hereafter Paper II). In the following subsections, we 
present the demonstration of the corresponding new formulas. This section is a bit technical and the confident reader 
can skip the demonstrations and go directly to subsection [l.6| . 

4-.1. A first simple approach 

We begin with a very simple model, in which the galactic disc is an infinitely thin disc extending to infinity, embedded 
in an infinite diffusive volume. These assumptions are not unreasonable as propagation of radioactive species is a local 
phenomenon and is not expected to be much affected by boundaries. We also ignore destructive spallations, galactic 
wind and energy losses. At the origin of the galactic plane stands a hole of radius rhoie- 
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Consider first the nuclei coming from spallations at time to and at a source point tq- They act as an instantaneous 
source term fnS{r — ro)^(t — to). The number that would reach the center at time t > to, with no decay, satisfies 



dt-^''f=^d-r['i-r)=f^'^''-'^^'^'-'^^ 
A solution of this equation is given by 

f{r -ro,t- to) = : exp \ - 



(MK{t-to))'^' I ^K{t-to). 
For radioactive species with a lifetime r = 7T0 (where tq is the rest frame lifetime), decay must be taken into account 
and we have instead, at r = 0, 

/(0,t) oc — exp|-— ii— 1 X cxp l-i^l (4) 



(t-to)'/^ ^K{t-to), 
This function is actually the propagator of this diffusion problem. In the permanent regime, the total quantity of these 
nuclei at the center of the hole is obtained by summing the solutions for a distribution of point sources acting at all 
instants from to = —00 to t (now). It leads to an integration over space and time 

in our case, n{ro) = for vq < Vhoie and n(ro) = no for tq > Vhoie, so that 



TV oc no l^^^^ 2.ro dro dto j^-^ exp x oxp {-^} 



X cxp 



t-tn 



Integration over ro is easily performed 

i-ooVF^ "^l ^K{t-to) 
Introducing x = \/{t — to)/T, 

Jo I ^Kt 

This integral is given by 

r dxe^v{-^-x'] = ^e-^^ (5) 
Jo X J 2 

so that finally 

The central density due to the full disc would be given by rhoie = 0, so that, introducing Z^ad = VK^yro, 
Nrnoi. (r = 0) = A(,,„,^=o) (r = 0) X e"'^'-'"^/'-'' (6) 
We emphasize that in the problem treated here, the time-dependent resolution by means of propagators that we used 
in this section is equivalent to solving the corresponding stationary equation {df/dt = 0). This latter approach is 
more convenient for implementing the physical effects we have neglected here. It will be used in the next sections and 
as expected, the same behavior will be recovered. 



4-2. Analytical solutions without energy losses and reacceleration 

We now turn to a more realistic model, namely a cylindrical diffusion box of radius R and half-height L, and a matter 
disc of half-height h, considered as infinitely thin. The local subdensity has to be represented by a circular hole of 
radius rhoie located at the position of the Sun. Compared to the previous approach, it has the advantage of correctly 
taking into account the presence of boundaries. It also takes into account destructive spallations and galactic wind. It 
is actually the model used in Paper I and Paper II, with an additional hole of radius rhoie- 
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To make the problem tractable analytically, we make some simplifying hypotheses. First, to preserve the cylindrical 
symmetry, we consider that the Sun (and the hole) is located at the center of the disc. Second, we suppose that the 
density of the parent species is uniform all over the galactic disc, i.e. it does not depend on the spatial coordinate r. 
The validity of these assumptions relies on the fact that the propagation of radioactive species is a local process. We 
will discuss further these points in section 4.4. 



Diffusive volume 




Galactic disc 



Fig. 3. Schematic representation of the model. Left picture is the actual geometry of the problem, with a hole of 
radius rhoie in the disc, at the position of the Sun. The right picture represents the simplified geometry (cylindrical 
symmetry is preserved) which is shown to give the same results. 



We want to solve the diffusion equation for a radioactive secondary (no primary sources) 

KAN''"'(r,z)~Vc V ; ; -rl^N""'(r,z) + 2hS(z) [rP™''(r)iVP'"'(r, 0) - F™"' (r)iV""''(r, 0)1 =0 (7) 

oz 

To keep notations simple, we consider only one parent nucleus. It would be straightforward to generalize to several 
parents. Here, iV""'* denotes the unstable nucleus, N^""^ the parent nucleus and F'"'"'^, F™*^' and denote respectively 
the production rate (by spallation) of the parent N^""^ into the radioactive nucleus iV""** , the destruction rate of iV""'* 
by inelastic scattering and its /3-disintegration rate. This equation differs from the no-hole case (equation (Al) of 
Paper I) because the terms F^"'°'^(r) — nLjsM{r)vaP^°'^ and F™'^'(r) — nLisM{f')vcr™'^^ now depend explicitly on r via 
the local interstellar density which reads 

riLiSMir) = 0(r - rhoie)nisM (8) 
where Q is the Heavyside distribution. Thus, we can rewrite equation in the form 

KAN'^^^r, z)-Vc -A^-T^^N'^^^r, z)+2hS(z) fFP'-°''e(r - rhoie)NP'''- (r, 0) - F™'^'e(r - rhoie)N""'' (r, 0)1 = 0(9) 

oz '■ ■' 

As for the no- hole model (see Paper I), a solution is found by expanding the density over the Bessel functions Jo{Cir/R) 
where the Q are the zeros of Jq. The unknown quantities iV""'*(r, z) are expanded as 

7V™-(r, z) = Y. Nr^i^Uo (Q^) (10) 

■i 

The known quantities are also expanded as 

^Yl^t{z)JoiCi^) with T,{z) ^ pT{p,z)Jo{Qp)dp 
where p — r/R. 

We now have to differentiate the case of spallation from that of destruction. For the first one, we can assume that the 
density of the parent depends very smoothly on the coordinate r so that we can safely take NP°'^{r, 0) = N^""^ {Rq,0); 
this approximation will be checked in ( |4.4| ). We must compute the Bessel transform of the Heavyside function 

= ^®(^-^)'^°(^»^)^^ 

With the notation x = Qp and using the fact that a primitive of xJo{x) is xJi{x), it can be shown that 



e,, = 



Jite)-^Ji(c.^) (11) 
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For the destructive spallation term, we have to Bessel develop the function Q{r — rhoie)N'"^''ir, 0). At variance with 
the precedent case, the distribution iV""*(r, 0) is expected to vary significantly across the hole so that the previous 
approximation cannot be made. Wc thus have 

oo 

e(r - r;,„,e)iV™^(r, 0) = J] f^r^o (C.^ 
where 



i=l 



piv"-(p,o)e(p-^ 



■MCiP) dp 



Inserting the Bessel development of N""-^ (r, z — 0) into this relation (see equation nG) , we find 
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Using the property 

pMQp)Jo{Cip) dp 



C - C 



<iPMQp)Ji{CiP) + C,oPJi{Qp)Jo{<^rP)\ for i ^ j 



pJ^iQp) dp = -p' [j2(C,p) + jiiQp)] 



else 



it follows that 



' hole 



E 



iV™^(0) 



R 



drfu 



R 

Ji 



-J! 

C,iT}iole 

R 



lole 



R 



(12) 



R 



Jo 



R 



Finally, putting everything altogether, the Bessel transform of equation (H) reads 



_ Vc_d_ 

dz^ K dz 



R^ K ) 







(13) 



where 9^ and O""* are given by equations (|ll| ) and (jlj). This equation is now very similar to equation (A6) of Paper 
I, and resolution proceeds as exposed therein. The solution in the halo is 



Nr\z) oc e^=^/2^ X sinh l^i^^— ^1 



where Si and Ai were defined in section ( 3.1.2 ) . The solution in the disc is obtained by integrating equation ([^) across 
the disc. 



2m 



'(0)-2iVr^(0)-^- 



K 



-a' 



K 



-e,iv?'"'(i?o,o) = 



(14) 



which gives 

^-(0) = e. X 2/.r^-''iv;-(i^o,o) _ 2^ 

The first term is very similar to the secondary source contribution of a no-hole model, but with an additional Qi factor 
taking correctly into account the effect of the hole on production spallations. The second term takes into account the 
effect of the hole on destructive spallations. It is expressed as 



hole 



^^^"(0) 



drhole 



R 



2rhoU 



7V™"(0) 



E 



Cjrhole 

R 



Ji 



R 

R 



(15) 



R 



Jo 



Cirh, 



R 



As this term depends on A'J"*(0), ( [l4| ) is a coupled set of equations and may be tricky to solve in practice. 

Actually, as the effect of spallations is expected to be small, we use a perturbative method. We first compute the 
solution of equation (Q) without the last term of equation (^5|), 

ivr'<°^(o)=.e.x^^™'^''"^(^-o) 



A' 



with 

A' ^A, + 2hr 



inel ^hole 



R 



R 
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The term we have neglected can now be taken into account by substituting the Nj-"''^{0) in equation (|15|) by these 
zero-order solutions 7V™*'*-°\ 



N 



uns,{0) 



Ji{Ci)A^ R 



(0) 



I ^ I '^l 



Cj'^hole 1 J- 



Ci'^hole 

R 



The new quantities 7V™'*''-^^(0) can then be used instead of 7V™'*'*-°^ (0) to estimate the correction at the next order. 
This procedure is repeated until convergence is reached. In practice, convergence is very quick and we only need to 
iterate twice. 



4-3. Behavior of the solution as a function of the hole radius 

A numerical study of this solution (see figure ^ shows that the dependence in rhoie can be expressed as, with good 
accuracy. 



iV,(0) «exp 



with Irad = \JK{E) 



(16) 



Table 5. Values of Irad for several radioactive nuclei and energies. 





TO (Myr) 


Ek = 100 MeV/n 


1 GeV/nuc 


10 GeV/nuc 




2.17 


0.075 kpc 


0.22 kpc 


0.95 kpc 


26^1 


1.31 


0.037 kpc 


0.11 kpc 


0.47 kpc 


36 Q 


0.443 


0.019 kpc 


0.056 kpc 


0.25 kpc 



The physical meaning of this dependence is better seen with the first approach (section 4.1). The relevant quantities 
to compare are the size of the hole and Irad, which represents the typical distance travelled by a radioactive nucleus 
before it decays. This distance increases with energy, because of time dilation and because high energy nuclei diffuse 
more efficiently (K grows with energy) . For a diffusion coefficient of the form K ~ KqPTZ^ , Irad can be expressed as 



^(5-l)/4 



■p-^k,nuc ) 



(17) 



Typical values are given at several energies and for several radioactive nuclei for Kq — 0.033 kpc^ Myr ^ and 6 = 0.6 
(a good model taken from Paper I) in Table |[ 



4-.4- Discussion of the validity of the hypotheses 

The most a priori questionable assumption is that we located the hole at the center of the galactic disc, though we 
know that the Sun is at a galactocentric distance of about R = 8 kpc. Figure ^ displays the radial distribution of a 
radioactive species in the disc for several energies. We see this distribution is only perturbed locally by the presence 
of the hole: at a distance of a few Irad, the flux is no longer affected by the hole. Conversely, the flux at the center 
of the hole is not affected by the boundaries of the diffusive box if they are farther than a few Irad- We could even 
consider that the galactic disc is an infinite plane. This means that the same distribution - and the same reduction 
factor exp{—r hole /Irad) - is obtained whatever the position of the hole, as long as it is far enough from the edges of 
the box (compared to Irad), which is the case for us. The density of parent nuclei must obviously be estimated at the 
position of the Sun (and not at the center of the galaxy !). 

We also assumed that the parent has a homogeneous distribution all over the disc, which is known to be wrong: 
the parent fiux is greater around the position of the sources (i? ~ 4 kpc, see for example Paper I, section 3.5.1 and 
references therein) and decreases regularly to zero at the edge of the disc {R = 20 kpc). However, as the radioactive 
propagation is a local phenomenon, at the scale of Irad, it is a good approximation as long as the parent distribution 
does not vary much on the spatial scale of Irad- Actually, even if the parent distribution varies linearly with radius, 
the approximation remains good. The probability dV that a radioactive nucleus detected at the center of the hole 
has been created by a spallation at a distance between r and r + dr can be computed explicitly in the simple model 



described in subsection 4.1, and it is found that (see Appendix for the demonstration) 



c?7'(emitted between r and r + drjdetected at center) oc 6(r — r/io/e) A^parent(^) e ''/''■"■^ dr (18) 
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Fig. 4. The radioactive flux is computed for ^® Al and ^°Be, for several energies (corresponding to several values of 
lrad{E)) and for several hole radii r^oie- This plot displays each flux divided by the flux obtained in the homogeneous 
case {rhoie = 0)1 a, function of rhoie/lrad- The solid line represents the function exTp{—rhoie/lrad)- 



so that most of the radioactive nuclei come from a ring extending from the edge of the hole (r = rhoie) to a few Irad 
away. As a result, if the parent nucleus density is not uniform, the spallation rate is determined by the effective density 
obtained as the average over the disc with the appropriate weight, given in the expression A.l (see Appendix) 

j^effecUve ^ II ^'^^ N^grentjr) Q{r - rhole) exp{-r/lrad) 
parent JJ N^^rentir) e(r - r^^i,) eM-r/lrad) 

We have compared this effective parent density to the actual parent density Nparent{RQ, z = 0) at the location of the 
Sun. For hole radii rtoie < 400 pc and for Irad < 5 kpc, the difference is less than one percent at 10 GeV/nuc, and is 
completely negligible for E <1 GeV/nuc. 



4-5. Energy losses and reacceleration 



In homogeneous models, the energy changes induced by energy losses and reacceleration are described by the equation 
(see Paper I and Paper II) 



AiN, + 2h 



dE \ 



b{E)N,{E) 



m 

■ dE 



(0) 



(19) 



where Ni is the Bessel transform of the final flux and iVj^"-* is the initial flux (before energy losses and reacceleration 
are applied). Because of the hole, the energy loss term b{E) now has a spatial dependence. It has the same value as 
above, except in the hole where there is no matter to interact with, so that 
bir,E) = Q{r-rhoie)biE) 

Thus, the Bessel transform of the quantity 8(r — rfioie)N{r, E), denoted I5i, must be introduced. As the spatial 
dependence of this function is exactly the same as the spallation term of the previous section, it is straightforward to 
write (see eq. n2) 



ME) 



N^{E) 



1 



hole 



2f/-.\ (-2 



R 

Qrho 
R 



Ji 



ole 



R 

Ci^hole 



(20) 



R 



+ OJi 



R 



Jq 



Ci'^hole 

R 



so that b{E)Ni(E) in ( p^ has to be replaced by b{E)l3i{E). It is again a coupled set of equations. The same perturbative 
approach as before is used: the new quantities N^"^'' affected by energy losses are estimated by replacing Ni by Nj^^\ 
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Fig. 5. Radial distribution of a radioactive species in the disc, across the hole, for rhoie = 200 pc. Numerical values 
were taken for ^^Al (radioactive) and ^^Si (stable). The distribution of radioactive is very sensitive to the presence of 
the hole. However this effect is local and vanishes at several rhoie- On the contrary, the distribution of stable species 
is not much affected. 

The new set N^^^ is then inserted instead of N^^\ etc. . .until convergence is reached. It is found that even in the 
presence of energy losses and reacceleration, the dependence on the hole radius is still very well described by (|l6|) . 

We note that reacceleration has a negligible effect on radioactive propagation, which may be understood as it takes 
much longer than the lifetime jtq to significantly reaccelerate a nucleus. 

4-. 6. Conclusion and inclusion in the propagation code 

The presence of a hole of radius r/ioie lowers the radioactive flux by a factor given with a good precision by 
ex.p{-rhoie/lrad) wheie Irad is given in (|l^). 

It is not possible in practice to directly use the analytical formulae described above because the double sums over 
Bessel indices are far too time consuming. Indeed, the accurate description of the hole requires a development of all 
functions over at least the 3> R/vhoieT^ ~ 100 first functions Jo{(ix). We used N = 5000, for which a very good 
convergence of the Bessel series was achieved, and the double summation takes N'^ = 2.5 x 10'' elemental steps for 
each iteration and each model. The computation time required to apply these formulas to all the models would have 
been too large. 

Thus, we preferred to take advantage of the exponential dependence discussed above, and the flux of a given radioac- 
tive nucleus in the presence of a hole is obtained from each model with no hole by a multiplication by exTp{— rhoie / Irad) ■ 

5. Experimental data and configurations of parameter space 

We now turn to the more standard diffusion aspects of propagation. The main ingredients of our diffusion model have 
been widely depicted in Paper I, and we will not be exhaustive about them here. Instead, we simply motivate the 
choice of configurations that are used in this paper. We also review the experimental data that we compared to our 
calculations. 

5.1. Sets of configurations used for the analysis 

Propagation is assumed to be a diffusive process occurring both in the galactic disc and in a halo of thickness L. This 
process is characterized by an energy-dependent diffusion coefficient of the form K{E) ~ Kq(3TZ^ where and 5 
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are parameters of the model and TZ is the rigidity. Propagation is affected by a galactic wind Vc perpendicular to the 
galactic plane and by the presence of Alfven waves of velocity Va- These are the five parameters of the model. All the 
sets of parameters consistent with measurements of stable nuclei 

analysis have been extracted and discussed in Paper I (see this paper for a detailed presentation of the model and 
more specifically figures 7 and 8). These sets of parameters are also consistent with the antiproton spectrum, as shown 
in Paper II. 

Some further cuts in our initial sets of parameters could probably be applied on physical grounds (see for example 
the discussion on the Alfvenic velocity Va in section 5.2, Paper II). We adopted a conservative attitude and used the 
whole set, as the aim of this paper is precisely to explore whether it is possible to obtain further constraints on these 
parameters from the study of radioactive species. 

5.2. The data set we used 

Several experiments in the last twenty or thirty years have measured radioactive isotopes in cosmic rays with increasing 
precision, at energies of a few hundred MeV/nuc. The early data - usually presented as the ratio of the radioactive 
isotope to its stable companion(s) - were affected by errors of around 25-30%. The latest published data have error 
bars reduced by a factor of two or three. In the following, we implicitly refer to three satellite experiments, namely 
Voyager, Ulysses and ace. Other experiments will sometimes be shown on figures but they will be purely illustrative 
since their accuracy is far smaller. 

The best measured ratio is probably ^"Be/^Be which corresponds to the lowest Z /3-radioactive nucleus. Data from 
Ulysses (Connell, 1998) and from ACE (Binns et al., 1999) are consistent, the quoted error bars being smaller for ACE. 
They are also consistent with the Voyager data point (Lukasiak et al., 1999) for which the quoted error is larger. We 
do not use the SMILI data, as the possibility that they are due to a statistical fiuctuation is not ruled out (Ahlen et 
al., 2000). 

As regards the radioactive chlorine isotope ^^Cl, results are usually provided as ^^Cl to total CI ratio. The only 
available data, to our knowledge, are those from Ulysses (Connell et al., 1998) with a la error of about 35%, and ACE 
(Binns et al., 1999) whose errors (even taken at Scr) are completely included in the Ulysses l-tr upper error band. 

We finally end with the ^^Al/^^Al ratio. Contrary to other radioactive ratios, the measurements show more prob- 
lems. Indeed, la data from Ulysses (Simpson and Connell, 1998) and ACE (Binns et al., 1999) exclude each other (the 
ACE central point is much lower than Ulysses'). Enlarging ACE error bars (which are smaller than Ulysses) to Str does 
not improve the compatibility. On the other hand, Ulysses data are fully compatible with la Voyager data (Lukasiak 
et al., 1994), whose uncertainty is still much greater than for the other two experiments. The possible discrepancy 
between some of these data will be addressed later. 

6. Results 

We now present our analysis of the radioactive nuclei abundances. First, in the diffusion model with no hole (hereafter 
the homogeneous model), we can compute the expected flux of radioactive nuclei for each set of diffusion parameters 
given by the analysis of stable nuclei (Paper I) and compare to the data. Then, the presence of a local underdense region 
discussed above is tested, and the radius rhoie of this region is introduced as an additional parameter. The corresponding 
models will be called inhomogeneous models. As the stable nuclei are almost not affected by the presence of the hole 
(see figure the sets of diffusion parameters given in Paper I are also used in the inhomogeneous case. 

To sum up the procedure, the sets of diffusion parameters {Kq, L, 5, Va and Vc) given in Paper I are used to 
compute the radioactive nuclei fluxes, for different hole radii rhoie- The special case rhoie = corresponds to the 
homogeneous models. 

We proceed as follows. We flrst focus independently on the ^"Be/^Be and ■^^Cl/Cl ratios. They are computed for all 
the sets of diffusion parameters compatible with B/C (given in Paper I) and for several hole radii rhoie - The constraints 
they induce on the parameters are studied. We then analyze simultaneously the two ratios ^''Be/^Be and ^^Cl/Cl and 
finally the three ratios ^OBe/^Be, ^eci/Cl and ^^k\/^'^k\. 

6.1. Constraints from Be Be 
6.1.1. Comparison to ACE data 

To begin with, we consider only the ^°Be/^Be ratio, for which we have three compatible measurements. This is also 
the ratio for which an accurate spectrum is likely to be available in a near future (by ISOMAX, AMS, . . . ). This ratio 
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is computed for each of the 5-parameter sets compatible with B/C given in Paper I, and for several hole radii rhoie- 
The result is compared to the ACE data and the parameters giving a ratio falling out of the 3cr ACE error bars are 
discarded. We also do the same with a more stringent limit of Icr, and the results are shown in the left panels of 
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Fig. 6. Representation of the models compatible with combination of B/C (see Paper I) and various radioactive ACE 
data. See text and notations in the figure. Upper panels display homogeneous models [rtoie = 0) in the plane L-5. 
Lower panels display inhomogeneous models [vhoie > 0) in the plane rhoie~5. 



figure |g, where the scatter of the configurations is plotted in the plane S-L for the homogenous case, and 5-rhoie in 
the inhomogeneous one. For illustrative purpose, we also show on the right panels of figure ^ the result of the same 
procedure with 3cr and la •^^Cl/Cl ACE data. The dots in this figure are merely the models obtained in the Paper I 
analysis. 

The upper figure shows that for homogeneous models, the ^''Be/^Be ratio alone further constrains L to a smaller 
region of the parameter space. For inhomogeneous models (lower part), the hole radius is constrained to values 
Thole ^ 100 pc. This can be easily understood: the value of ^''Be/^Be in an inhomogeneous model is given by the 
corresponding value in a homogeneous model, decreased by the exponential factor exp(— r/ioie/^rad). As the initial 
homogeneous parameter sets compatible with B/C give a wide range of ^"Be/^Be values, those which are larger 
than the data are redeemed in inhomogeneous models. For holes larger than 100 pc, the exponential decrease is too 
important and theoretical predictions are too low to fit ACE data. 

If we do not believe in the presence of the bubble, then our conclusions are similar to previous works: higher values 
of L are preferred in homogeneous models {L'>^ A kpc). For inhomogeneous models (including ruoie — 0), the allowed 
ranges for the other diffusion parameters are not much affected if compared to the results of the stable nuclei analysis 
(Paper I). The ■^^Cl/Cl ratio yields similar conclusions. 
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6.1.2. Discussion 

Here we want to discuss qualitatively two other important points: the spectral behaviour and the interrelation of the 
three radioactive ratios. It will allow a more intuitive grasp of the combined analysis presented in the next section. 

We first select all the homogeneous models compatible with ^°Be/^Be {la and 3cr ACE data). Then, for each point 
of the resulting subset, the full propagation code is run and for the three radioactive species, maximal and minimal 
ratios are calculated at all energies and generate an envelope. These envelopes are represented in the upper part of 
figure ^ They contain all the allowed values of the above-mentioned ratios which are "compatible with B/C" and 
"compatible with ^°Be/^Be ACE data la" . The same procedure is applied to the inhomogeneous models rhoie = 80 pc 
(this value will be favored in the combined analysis). 




Fig. 7. Envelopes of the spectra obtained with all the models compatible with ^°Be/^Be ACE la (upper) and 3a 
(lower) for the three radioactive species. Solid lines are for homogeneous models {vhoie — pc) and dashed lines are 
for inhomogeneous models [rhoie — 80 pc). Data in the lower panel are from ACE (circles), Ulysses (crosses). Voyager 
(filled squares) and isee (filled triangles, Wiedenbeck, 1985) 



We now want to address the following questions: (i) is it possible, from the radioactive spectra, to experimentally 
distinguish inhomogeneous from homogeneous models, (ii) is it possible to break the degeneracy in the propagation 
parameters, and (iii) what would be a clear signature of the presence of the hole ? 

As regards the first question, the upper left panel (^"Be/^Be) clearly shows that the answer is no, as even for 
extreme models {ruoie ~ 80 pc), there is still an overlap between the homogeneous and inhomogeneous envelopes. 
Concerning question (ii), a partial answer has been given in the previous section: the ^"Be/^Be ratio is degenerate in 
the diffusion parameters at the ACE energy. Figure ^ shows that it is actually true for all energies. We finally turn to 
question (iii). From left to right panels of figure |^, the nuclei have smaller and smaller lifetimes, so that the influence 
of the hole is larger and larger. In particular, the ratio ^^Cl/Cl is the most sensitive to rhoie, and at low energy, 
homogeneous and inhomogeneous envelopes (for rhoie = 80 pc) do not overlap. We thus expect that the combination 
of these ratio may break the degeneracy in rhoie- 
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We also show, in the lower part of figure |j, the envelopes obtained as above but with the weaker constraint that 
^°Be/®Be is in the Scr ACE error bars. For rhoie — 80 pc and at the ACE data energy, the ^°Be/^Be is never larger than 
0.12, which corresponds to the largest ^°Be/^Be obtained in the homogeneous models, decreased by the corresponding 
exponential factor. This is not shown in the figures, but for Vhoie = 50 pc, the lower panel dotted lines would be almost 
superimposed on the solid lines, whereas for rhoie — 100 pc, they would only delimitate a very narrow strip. We also 
see that inhomogeneous models {rhoie = 80 pc) are favored by the '^^Cl/Cl data. On the other hand, all the models 
compatible with ^"Be/^Be clearly fail to reproduce the 3cr ACE ^^A1/^''A1 data. 

6.2. Combined analysis of^"Be/^Be and '^^ CI/ CI 

We now analyze simultaneously the two ratios ^'^Be/^Be and ■^^Cl/Cl. From the initial set of parameters compatible 
with B/C, we select those giving values of ^°Be/^Be and ^^Cl/Cl in the error bars of ACE (Icr and 3cr). As expected 
from figure ^, there is no homogeneous model which fulfills the above condition at Icr (the two regions delimited by 
the filled circles in the left and right panels do not overlap). This can be seen in figure ^ which displays the models 
compatible with ACE at the Icr (filled circles) and 3cr (empty circles) level. We saw in the previous section that high 
values of L are independently favored by the two ratios ^"Be/^Be and ^^Cl/Cl for homogeneous models, and it is 
natural to recover this trend in the combined analysis, which points towards i ^ 6 kpc at the Scr level. The allowed 
range for the other parameters are not much changed. 

As regards the influence of rhoie, we see that the combined analysis naturally favors hole radii 60 pc < Vhoie ^ 80 pc 
(as seen before, values rhoie ^ 100 pc are independently excluded at the 3cr level), which is in full agreement with LISM 
observations. Note that for these particular models, values L > 12 kpc are excluded at the Icr level. 

6.3. Combined analysis all three radioactive: problem with data ? 

The next logical step is to include the ^^Al/^^Al ratio in the analysis. We first display in figure ^ the values of this 
ratio for the models discussed in the previous section, i.e. compatible with ^"Be/^Be and ^^Cl/Cl ACE data at 3cr. We 
also show the ^^Al/'^^Al experimental bounds from ACE and Ulysses at 3cr. 

The first strong conclusion is that all these models (homogeneous and inhomogeneous) give similar ^^Al/^^Al values 
at ^ 200 MeV/nuc. Moreover, it is seen that these values are not compatible with ACE data, even at the 6cr level. We 
will come back to this point later on. As pointed out in section other experiments are not compatible with the 
ACE data for this particular ratio. To consider the possibility that the ^^Al/^^Al value may be higher than hinted by 
the ACE measurement; we repeat the previous analysis applying Ulysses data for this ratio. 
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Fig. 8. Representation of the models compatible with B/C plus both -'^"Be/^Be and ■^^Cl/Cl ACE 3ct (open circles) 
and Icr (filled circles). Left panel displays homogeneous models {rhoie = 0) in the plane L-S. Right panel displays 
inhomogeneous models {rhoie > 0) in the plane rhoie~5- 
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Fig. 9. Representation of the ^^Al/^^Al ratio at the ACE energy as a function of L for the models compatible with 
(i) B/C in homogeneous models (small squares), (ii) B/C plus ^°Be/^Be and '^^Cl/Cl ACE 3cr (medium squares are 
for homogeneous models and big squares are for inhomogeneous models) . Solid lines represent the 3cr error band from 
ACE, dashed lines represent the Scr error band from Ulysses. 
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Fig. 10. Representation of the models compatible with B/C plus ^°Be/^Be and ^'^Cl/Cl ACE 3a plus ^''Al/^'^'Al Ulysses 
3ct (open circles) and la (filled circles) data. Left panel displays homogeneous models {rhoie = 0) in the plane L-S. 
Right panel displays inhomogeneous models {rhoie ^ 0) in the plane rhoie~S. 
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In figure [Tc|, we show the models compatible with the ^"Be/^Be and ^^Cl/Cl ACE 3a data, plus the ^^Al/^'^Al Ulysses 
3(7 data (open circles) and la data (filled circles). In the latter case, only inhomogeneous models with rhoie ~ 100 pc 
are consistent with the experimental values for the three ratios. In particular, we emphasize that homogeneous models 
are excluded. 
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Fig. 11. Envelopes of the spectra obtained with all the models compatible with ^°Be/^Be ACE and "^^Cl/Cl ACE 3a plus 
^^A1/^''A1 Ulysses 3a for the three radioactive species. Solid lines are for homogeneous models {rhoie = 0) and dashed 
lines are for inhomogeneous models {rhoie > 0). Data in the lower panel arc from ACE (circles), Ulysses (crosses). 



However, it must be stressed that for these "good" models, the ^°Be/^Be and "^^Cl/Cl ratios have the lowest 
possible value still compatible with 3a ACE data, whereas the ^^Al/^^Al ratio has the maximal value compatible with 
la Ulysses data, so that these ratio are only marginally consistent. In figure |ll|, we also present the envelopes obtained 
as in the previous section, but with the more stringent condition that the ^"Be/^Be and '^^Cl/Cl ratios are within the 
1 a ACE error bars. 

We did not make the analysis with the Voyager data point instead of Ulysses because the error bars are much 
larger. However, we notice that the central value provided by this experiment is in better agreement with our derived 
^^Al/^^Al ratios (see figure @). 



7. Conclusion and discussion 

We have presented the analysis of the three radioactive ratios -'^"Be/^Be, "^^Cl/Cl and ^^Al/^^Al in the framework of 
a diffusion model, taking into account the presence of a local underdensity of radius rhoie- We find that the ^"Be/^Be 
and ^^Cl/Cl ratios are compatible with ace data if rhoie ^ 60 — 80 pc, but that no model is simultaneously compatible 
with the ACE data for the three ratios. However, if we consider other ^^Al/^^Al data (Ulysses or Voyager), some 
models are marginally consistent and the cases 60 pc < rhoie ^ 100 pc are preferred. These values are in agreement 
with independent estimations of lism studies. The presence of this hole would also lead to an attenuation of the 
expected "'^^C flux of the order of 10"'''. Thus, it seems a very difficult task to detect any ^''C in cosmic radiation, 
unless some local source is present. For example, radioactive nuclei could be produced in the very local fluff (with 



a very low density, see section 3.2.1 ). This has been more realistically modelled by Ptuskin and Soutoul (1990) and 
Ptuskin et al. (1997) with a three-layer bubble. It would be straightforward to adapt our model to take into account 
a multi-shell contribution (as long as the symmetry is preserved). Moreover, all these models assume that the shells 
have the same center, which is surely wrong, but the generic behaviour in exp{— rhoie /Irad) is expected to hold even 
in a more complicated geometry. 

We have not yet discussed solar modulation. All the data have slightly different solar modulation parameters. We 
used the simple force field scheme to modulate our fluxes (as in Paper I and Paper II). We checked that even taking 
extreme values for the modulation parameters yields about 5% changes in calculated ratio. This does not affect the 
results of our analysis. 
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Our results point out an important feature of cosmic rays. Either there is a problem with the Al data, or there 
is a fundamental problem with the diffusive approach of radioactive nuclei propagation. Alternatively, there could 
be a problem with cross sections. Other similar studies (see for example Strong and Moskalenko, 1998 and further 
developments) should be able to confirm our results when incorporating a hole in their model. 
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Appendix A: demonstration of formula (^) 

In this section, we compute the probability dV that a radioactive nucleus detected at the center of the hole has been 



created by a spallation at a distance between r and r + dr, in the framework of the simple modelling of subsection 4.1 

It is a conditional probability, and as such it can be expressed, using Bayes theorem, by 

7^ r 1 I , 1-^1 ^ r , 1-^1 1 ^ dV {created at r} 

dV {created at r|detected m 0} = V {detected m 0|created at rj x ^ ^ — ^ ^ . — — 

cx V {detected in 0| created at r} x dV {created at r} 

The first term is given by the expression (0) 

r°° 1 t ^ dt 
V {detected in Olcreatcd at rj cx / cxp < — -— \ ^tttt 

and the second term is simply given by 







JT. r +^ + 1 Qir~rhoie)'2TrrN{r)dr at/\^ 

dV {created at rj = jr- — cx e(r - rhoie)rN{r) dr 

j 2TTrN[r) dr 

where N{r) is the density of parent nuclei. This gives 

/■°° ( r^ t ] dt 

dV {created at rldctected in 0} cx 9(r — rhoie) r N{r) dr x exp < > -^r-pr 

Jo { '^Kt 7ro J t-^i^ 

Defining x — r/\/AKt, the integral is proportional to 

exp < — X > dx 

^ \ AKyrox^ J 

The value of this integral is given by (^ which gives 

dp {created at r|detected in 0} cx &{r - r^oie) N(r) e"''/'''"'dr (A.l) 

The proportionality coefficient of this relation would easily be obtained by imposing that J^^p dV — 1. 
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